Engineering of oscillatory quantum states by parametric excitation 



o : 

<N ■ 



^ : 

Oh- 
^— > ! 

c : 

=5 

cr 



> 
in 

00 

cn 

o 



X 



T. V. Gevorgyan^B and G. Yu. Kryuchkyan 1 - 2 '!!] 

1 Institute for Physical Researches, National Academy of Sciences, 
Ashtarak-2, 0203, Ashtarak, Armenia 
2 Yerevan State University, Alex Manoogian 1, 0025, Yerevan, Armenia 

We consider preparation of nonclassical oscillatory states in a degenerate parametric oscillator 
combined with phase modulation. In this scheme intracavity oscillatory mode is excited by train of 
Gaussian laser pulses through degenerate down-conversion process and phase modulation element 
inserted in a cavity leads to anharmonicity of oscillatory mode. We demonstrate production of 
nonclassical oscillatory states with two-fold symmetry in phase-space including the superposition of 
Fock states and quantum localized states on the level of few excitation numbers and in over-transient 
dissipative regime. 
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I. INTRODUCTION 

There is currently a wide effort to construct various 
artificial nonlinear oscillatory systems showing quantum 
behavior. Such systems operated in quantum regime 
have become more important in both fundamental and 
applied sciences, particularly, for implementation of ba- 
sic quantum optical systems, in engineering of nonclas- 
sical states and quantum logic. In these systems the ef- 
ficiency of quantum effects requires a high nonlinearity 
with respect to dissipation. The nonlinearity breaks the 
equidistance of oscillatory energy levels that allows se- 
lectively excites the oscillatory states by monochromatic 
driving analogous to that take place for electronic states 
of atomic systems. Thus, the limited two-level and three- 
level systems can be realized approximately for oscilla- 
tory systems. 

An oscillator operated in quantum regime is naturally 
described by Fock states that are states with definite 
numbers of energy quanta. The preparation and use 
of Fock states and various superpositions of Fock states 
form the basis of quantum computation and communi- 
cations pj. However, excitations of oscillatory systems 
usually lead to production of coherent states nearly in- 
distinguishable from a classical state, but not quantum 
Fock states. In this reason quantum oscillatory states 
are usually prepared and manipulated by coupling oscil- 
lators to atomic systems. In this way, a classical pulse 
applied to the atomic states creates a quantum state that 
can subsequently be transferred to the harmonic oscilla- 
tor excited in a coherent state. The systematic procedure 
has been proposed in [2] and has been demonstrated for 
deterministic preparation of mechanical oscillatory Fock 
states with trapped ions Q and in cavity QED with Ry- 
dberg atoms [J] . Most recently the analogous procedure 
has been applied in solid-state circuit QED for determin- 
istic preparation of photon number states in a resonator 
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by interposing a highly nonlinear Josephson phase qubit 
between a superconductive resonator [5|. 

Recently, it has been shown that production of Fock 
states and Fock states superpositions or qubits can also 
be realized in over-transient regime of an anharmonic dis- 
sipative oscillator without any interactions with atomic 
and spin- 1/2 systems |6|. Preparation of the lower Fock 
states |1) and superposition state \ty) = ^(|0) + |1)) has 

been demonstrated for the lowest excitation |0) — > |1) 
with complete consideration of decoherence effects. For 
this goal the strong Kerr nonlinearity as well as the ex- 
citation of resolved lower oscillatory energy levels with 
specific train of Gaussian pulses have been considered. 

In continuation of this paper [6| here we propose the 
other oscillatory scheme based on two-quantum resolved 
excitations of oscillatory levels |n) — s- \n + 2). Thus, 
this scheme is the Duffing oscillator for the mode cre- 
ated in the process of degenerate down-conversion under 
an external field. One of the possible realizations of this 
scheme is the parametrically driven anharmonic oscilla- 
tor (PDAO) that consists from intracavity degenerate 
parametric oscillator combined with phase modulation 
element inserted in a cavity. In this scheme intracav- 
ity oscillatory mode is excited by strong field through 
degenerate down-conversion and phase modulation leads 
to anharmonicity of oscillatory mode. 

Note, that parametrically driven anharmonic oscillator 
based on cascaded parametric oscillator and third-order 
phase modulation has been proposed and studied in the 
papers 0] , [1] for the special cases without consideration 
of effects of dissipation and quantum fluctuations. Quan- 
tum theory of monochromatically driven parametric os- 
cillator combined by phase modulation has been devel- 
oped in terms of the Fokker-Planck equation in complex 
P representation Q, The exactly solution of this 

equation (that means consideration of all order of dissi- 
pation) has been obtained in steady-state regime and the 
Wigner function of oscillatory mode has been obtained 
in analytical form using this solution [9]. These results 
strongly demonstrate the vanishing of Fock states and 
quantum superposition states in over transient dissipa- 



2 



tive operational regime of PDAO within framework of 
the Wigner function that visualizes quantum effects as 
negative values in phase-space. Really, analytical results 
for the Wigner functions of PDAO under monochromatic 
excitation [9| has been calculated as positive in all ranges 
of phase space in the steady-state regime. 

In this paper, we demonstrate that in the specific 
pulsed regime of PDAO the production of nonclassical 
oscillatory states with two-fold symmetry in phase-space 
consisting form the Fock states |0), |2), as well as super- 
position of Fock states |\&) = 775 (|0) + |2)) are realized. 
These states can be created for time intervals exceeding 
the characteristic time of decoherence. 

The quantum regime of PDAO requires a high third- 
order, Kerr nonlinearity with respect to dissipation. The 
largest Kerr nonlinearities for oscillatory systems were 
proposed for cooling nano-electromechanical systems and 
nano-opto-mechanical systems based on various oscilla- 
tors [lj] , [H| • Superconducting devices based on the non- 
linearity of the Josephson junction (J J) that exhibits a 
wide variety quantum phenomena (l6j-[22l| offer an un- 
precedented high level of nonlinearity and low quantum 
noise. In some of these devices dynamics are analogous to 
those of a quantum particle in an oscillatory anharmonic 
potential 23 1, [24|]. Note, that comparison of third-order 
nonlinearities taking place for various quantum devices 
have been recently analyzed in [25| . 

The paper is arranged as follows. In Sec. II we describe 
parametrically driven nonlinear oscillator under pulsed 
excitation and descibe phase-space symmetry properties 
of the model. In Sec. Ill we shortly discuss PDAO under 
a monochromatic driving for both transient and steady- 
state regimes. In Sec. IV we consider production of non- 
classical oscillatory states in the pulsed regime of PDAO. 
We summarize our results in Sec. V. 



II. MODEL DESCRIPTION 

In this section we give the theoretical description of 
the system. The nonlinear oscillator driven parametri- 
cally by train of pulses and interacting with a reservoir 
is described by the following Hamiltonian 



H = HojQa + a + hx{a + a) 2 + 



m(E(t)e 



-iujt a +2 



+ E+{t)e iuJt a 2 ) + H lo 



(1) 



where a + , a are the oscillatory creation and annihilation 
operators, luq is the oscillatory frequency, x is the non- 
linearity strength proportional to the third-order suscep- 
tibility. The coupling constant Q is proportional to the 
second-order susceptibility and the time-dependent am- 
plitude of the driving field E(t) = Eof(t) consists from 
the Gaussian pulses with the duration T which are sep- 
arated by time intervals r 



Hi oss = aT + + a + T is responsible for the linear 
losses of oscillatory state, due to couplings with heat 
reservoir operators giving rise to the damping rate 
7. The reduced density operator p within the frame- 
work of the rotating-wave approximation, in the inter- 
action picture corresponding to the transformation p — > 
e -i{ui/2)a at p e i{u)/2)a at j g g 0vernec j \yy the master equa- 
tion 



dp 



-[H +H int ,p] 



E ( L ^ L t - \ L t^iP - IpLfL^) , (3) 



where L\ = \J (N + 1)70 and L2 — \JTFja + are the Lind- 
blad operators, 7 is a dissipation rate and N denotes the 
mean number of quanta of heath bath, 



Ho = HAa + a, 
H int = h X (a+a) 2 + hQ(E(t)a+ 2 + E(t)*a 2 ), (4) 

and A = luq — uj/2 is the detuning between half frequency 
of the driving field uj/2 and the oscillatory frequency uiq. 
Two last terms in the interaction Hamiltonian describe 
the self-phase modulation (SPM) of the oscillatory mode 
and the parametric three- wave interaction between semi- 
classical driving field and the oscillatory mode, respec- 
tively. 

To study the pure quantum effects we focus on the 
cases of very low reservoir temperatures which, however, 
ought to be still larger than the characteristic tempera- 
ture T 3> T cr — tvy/ke- This restriction implies that dis- 
sipative effects can be described self-consistently in the 
frame of the Linblad Eq. ([3]). For clarity, in our numer- 
ical calculation we choose the mean number of reservoir 
photons N = 0. 

It is evident that the system displays definite symme- 
try properties in phase-space. Really, considering the 
transformations 



H' = U^HU, p = U~ x pU 



with the unitary operator 

U = exp (ida + a) 



(5) 



(6) 



we verify that the interaction Hamiltonian §4$ satisfies 
the commutation relation 



[H, U] = 0. 



(7) 



The analogous symmetry takes place for the density op- 
erator of oscillatory mode 



/w = E< 



-(t-t„-nr) 2 /T 2 



(2) 



[p(t),U}=0. 



(8) 
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One of the most important conclusions of such symme- 
tries is related to the Wigner functions of the oscillatory 
mode 



1 



W (a) = -z / d^Tr pe 



(9) 



where we perform rotations by the angle 9 around the ori- 
gin in phase spaces of complex variables a corresponding 
to the field operators a in the positive P-representation. 
Indeed, in the polar coordinates r, 9 of the complex 
phase-space plane X — (a + a*)/2 — rcos9, Y 
n ft" jH = rsin# we derive that the Wigner func- 
tion displays two-fold symmetry in its rotation around 
the origin of phase-space 



W (r, 9 + tt) = W (r, 9) 



(10) 



It is well known that this symmetry take place for the 
degenerate optical parametric oscillator (OPO) and re- 
flects on the phase-locking phenomenon in above thresh- 
old regime of OPO. According to phase-locking the mode 
of sub-harmonic generated in OPO are produced with 
well-defined two phases 27]. The result (10) shows that 
this situation takes place also for the combined system 
under consideration. The illustrations will be presented 
below on the Figs. [4j [6] 

In the following the distribution of oscillatory excita- 
tion states P(n) = (n\p\n) as well as the Wigner func- 
tions 



w(r,e) 



.Pnm(t)W mn (r,9) 



(11) 



in terms of the matrix elements p nm — (n\p\m) of the 
density operator in the Fock state representation will 
be calculated. Here the coefficients W mn (r, 9) are the 
Fourier transform of matrix elements of the Wigner char- 
acteristic function. 

In this paper, we demonstrate that in the specific 
pulsed regime of PDAO and for the case of strong nonlin- 
earities the Fock states as well as superposition of Fock 
states are realized . These states can be created for time 
intervals exceeding the characteristic time of decoher- 
ence. The corresponding Wigner functions of oscillatory 
mode show ranges of negative values and gradually de- 
viate from the Wigner function of an PDAO driven by 
monochromatic driving. 

Nevertheless, the time evolution of PDAO driven by a 
coherent force cannot be solved analytically for arbitrary 
evolution times and suitable numerical methods have to 
be used. We solve the master equation Eq. ([3]) numeri- 
cally based on quantum state diffusion method [28]. The 
applications of this method for studies of the driven non- 
linear oscillators and the parametric optical oscillator can 
be found in [29|-[3lj. In the calculations a finite bases 
of number states \n) is kept large enough (with n max is 
typically 50) so that the highest energy states are never 
populated appreciably. 

It should be also noted that the investigation of quan- 
tum dynamics of a driven dissipative nonlinear oscillator 



for nonstationary cases, particularly, for various pulsed 
regimes, is much more complicated and only a few papers 
have been done in this field up to now. Quantum effects 
in nonlinear dissipative oscillator with time-modulated 
driving force have been studied in the series of the papers 
29]-[32fl in the context of a quantum stochastic resonance 
28 1 , quantum dissipative chaos (30j , l3~Tll , and quantum 
interference assisted by a bistability [32| . 



III. PDAO IN MONOCHROMATIC FIELD 

In this section, we shortly discuss the combined dis- 
sipative oscillator under monochromatic excitation, con- 
sidering f(t) = 1 in the Hamiltonian (f5]). At first, we 
present the results in the semiclassical approach and the 
standard linear stability analysis with respect to small 
deviations from steady state in terms of the amplitude of 



the oscillatory mode a 



,1/2 



exp 



(itp) [9j. The intensity 
n (in photon number units) and the phase of the mode 
ip in stable above-threshold regime is determined by the 
following expressions: 



n = f[^ + (J-l)n 

2x 7 

sin($ - 2tp) = J~ 1/2 , 



(12) 



where J 



-I and $ is the phase of the driving field 



Eq = I 1 ' 2 exp(i§). Above threshold regime takes place 
for I > I t h = qt(1 + t-) and the regular behaviour is re- 
alized for negative detunings. These results are obtained 
in over transient regime and for large oscillatory mean 
excitation numbers n >> 1. 

We now present results for strong quantum regime of 
PDNO that is realized for x/l > 1. In the absence of any 
driving, the quantized vibration states of nonlinear oscil- 
lator are the Fock states \n) which are spaced in energy 
E n = Eo + hioon + hxn 2 with n = 0, 1, .... The levels form 
a anharmonic ladder (see, Fig. HJa)) with anharmonicity 
that is given by E21 — E10 — 2h\- Below we concen- 
trate on quantum regimes for the parameters when os- 
cillatory energy levels are well resolved considering near 
to resonant two-photon transitions between lower num- 
ber states |0) — > |2), |1) — > |3). Thus, we assume non 
strong excitation regime Efl/^f = 7 and that the detuning 
^2 = -^(-E^ — -Eo)— w = 2A+4x meets the near to resonant 
condition, 5 2 = I67. The detuning between the transi- 
tion |1) |3) reads as S3 = \{E 3 -E x )-bJ = 2A+8x and 
for this parameters is more large than 62, i.e. £3 = 367. 
For comparison we add in Fig. [ljb) also the lower levels 
and the oscillatory transitions for the scheme considered 
in 

Now, we present on Fig. [5] the numerical results for 
monochromatically driven PDAO in dependence of the 
parameters: X, A, Efl. Transition the system to the 
steady state for time intervals t ^> 7" 1 is depicted in Figs. 
[5] (a,b) on the time-dependent mean excitation number 
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n=0 



(a) 



8 n=2 
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FIG. 1. Transitions between energetic levels of anharmonic 
oscillator: parametric driving with two-quanta detunings A = 
cj — w/2 and 82 = j(^2 — E a ) — uj — 2A + A\ (a); standard 
driving with one-quanta detunigs A = Wo — u) and <5i = 4 (_Ei — 
E )-w = A + x (b). 
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(c) 




FIG. 2. Time evolution of the averaged excitation numbers 
(a); the Rabi oscillations of the state populations with deco- 
herence which suppresses beating (b); the distribution of ex- 
citation numbers (c); and the Wigner function (d) for PDAO 
in steady-state regime. The parameters are: A/7 = —2, 
Xh = 5, MV7 = 7 




FIG. 3. Time-evolution of averaged excitation number (a) 
and populations of Fock states (b). Time-dependent struc- 
ture of pulses is indicated below n(t) in arbitrary units. The 
parameters are: A/7 = —2, xll = 5, Efl/y = 10, r = 47 -1 , 
T = 0.57" 1 



that the steady-state solution of the Fokker-Planck equa- 
tion has been found using the approximation method of 
potential equations [33| . The validity of this solution has 
not been checked in the strong quantum regime of PDAO 
operated on the level of few excitation number. On this 
reason we calculate the Wigner function numerically on 
the base of numerical simulation of master equation by 
using quantum state diffusion method [281 ] . The results 
are displayed in the Figs[2] (d). As we see, the Wigner 
function displays two humps and has two-fold symmetry 
in phase-space under the rotation on angle 7r around its 
origin (10). This effect reflects the well-known phenom- 
ena of phase locking that takes place for both OPO and 
OPO combined with phase modulation. In semiclassical 
approach according to the equation (12) phase-locking 
in above threshold regime means the forming of two sta- 
ble states of oscillatory mode with equal photon numbers, 
but with two different phases which are $/2 and $/2 + 7r. 
In quantum treatment of PDAO this phenomenon is dis- 
played as two-humped structure of the Wigner function 
as is has been demonstrated analytically in [lfj. Here 
we demonstrate that similar effect occurs also for PDAO 
in strong quantum regime on the level of few excitation 
numbers. 



IV. PRODUCTION OF OSCILLATORY 
QUANTUM STATES IN THE PULSED REGIME 



and the populations of Fock states. As we noted above, 
steady state results, in that number probability distri- 
bution of oscillatory excitation numbers and the Wigner 
function of oscillatory mode have been obtained analyti- 
cally in terms of the exact solution of the Fokker-Planck 
equation Q, [l(| • Particularly, the solution for the 
Wigner function of oscillatory mode involving quantum 
noise in all order of perturbation theory and in steady- 
state regime is positive in all phase space and hence does 
not describe Fock states in over-transient regime. Note, 



Now we are able to present the results concerning pro- 
duction of oscillatory quantum states for pulsed excita- 
tion regime of PDAO. Since in this paper we focus on a 
system that remains close to its ground, vacuum state, 
we assume a small number of excitation number in the os- 
cillator. We consider interaction time intervals exceeding 
the characteristic time of dissipative processes, t 3> 7 _1 , 
however, in the nonstationary regime that is conditioned 
by the specific form of pulsed excitation. If the paramet- 
ric excitation stipulated by train of the Gaussian pulses 
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the ensemble-averaged mean oscillatory excitation num- 
bers, the populations of oscillatory states and the Wigner 
functions are nonstationary and exhibit a periodic time 
dependent behavior, i.e. repeat the periodicity of the 
laser pulses in an over transient regime. 

As it has been shown in Q the most striking signature 
of periodically pulsed NDO is the appearance of various 
quantum states involving Fock states and interference 
between Fock states in over transient regime. Excita- 
tion due to parametric interaction leads to production of 
new class of oscillatory states that have not been possi- 
ble to generate in standard way by one-photon excita- 
tion. These states include squeezed states or superposi- 
tions of Fock state, as well as quantum localized states 
on the level of few excitation numbers. We demonstrate 
this point for near to two-quanta resonant transitions 
|n) — > |n+2) for the oscillatory parameters used above in 
Figs. [3] and S] however, for the periodically pulsed regime 
with the duration T — O.57 -1 , and the time interval 
between pulses t = Aj^ 1 . The evolution of the aver- 
aged excitation numbers and the population of the lowest 
states |0), |1) and |2) are depicted in Figs. [3]while time- 
dependent structure of pulses is indicated below n(t) in 
arbitrary units. As we see, for over transient regime the 
time-modulation of the averaged excitation number and 
the populations of oscillatory states repeat the periodic- 
ity of the pump laser. Besides these results, the excita- 
tion numbers for the definite time- intervals during pulses 
and the corresponding Wigner functions of oscillatory 
mode are shown in Figs. |3J As we see the evolution of an- 
harmonic oscillator under parametric excitation involves 
number of quantum oscillatory states which are visual- 
ized on the Wigner functions for typical five measurement 
time-intervals. We indicate the choosing of these time- 
intervals on the curve of main excitation number (see, 
FigEl (a)). For t = kr - 2T, k = 2, 3, ... the oscillatory 
mode is in the vacuum state, therefore the corresponding 
Wigner function is Gaussian (see, Fig[J](f)). After evo- 
lution of the system, for the time intervals t = kr — 1.8T, 
the state is squeezed in phase space. The typical result 
characterizing parametric double excitation of oscillatory 
mode is effective simultaneous production of |0) and |2) 
states that are depicted in Fig. [4] (c) for the time in- 
tervals t = kr — OAT. The populations of these states 
is approximately equal one to the other and the corre- 
sponding Wigner function (see, Fig0] (h)) demonstrates 
the interference fringes on phase-space between |0) and 
1 2) states. It is easy understand that this Wigner func- 
tion approximately describes the pure superposition state 
I Vfr) = _^(|0) + |2)), i.e. for these time-interval PDAO 

produces superposition of Fock states |0) and |2). It is im- 
portant that this quantum interference effect is realized 
when the oscillatory mean excitation number reaches its 
maximal value. For further increasing of time interval, 
i.e. for t — kr, the interference fringes on phase-space are 
deformed and the Wigner function consists from two lo- 
calized peaks with the ranges of negativity between them 
(see, FigfJ] (d, i)). Then, near to the end of the pulses, 



i.e. t = kr + 0.6T, the quantum interference is vanished 
and production of localized state takes place for PDAO 
(see, Fig|4](e,f)) as has been obtained in the steady state 
regime. It should be noted that these results are in accor- 
dance with the phase symmetry properties of PDAO (10). 
We conclude that production of quantum interference is 
realized in the vicinity of time-intervals where popula- 
tions of |0) and |2) states are crossed (see, Fig[3Jb)) in 
over transient regime. We assume that the control of de- 
coherence in this case take place due to application of 
suitable tailored, synchronized pulses (see, for example, 
[34j . (HI). Indeed, quantum interference is realized here 
if a mutual influence of pulses is essential, (for t/T = 3.14 
on Fig. EJ. 

The next regime that we consider is that where the 
pulsed excitation is tuned to the exact resonance, S — 
2 A + 4% = 0. It is evident that in this case the Fock 
state 1 2) can be effectively produced if a low excitation 
is used. Similarly to what was done above in dispersive 
regime, we consider at first the evolution of averaged ex- 
citation numbers as well as the population of the lowest 
states |0), |1) and |2) by two-quanta excitation. The 
results depicted in Figs, fallows us to choose the time- 
intervals within a pulse for which the maximal probabil- 
ity of production of |2) Fock state is realized. The time 
evolution of the probabilities Pq, Pi and P2 of |0), |1) and 
1 2) states starting from the vacuum state are depicted in 
Fig. OJb). As we see, that the maximal weight 0.6 of P2 
is realized for the definite time-intervals of measurement 
t = kr - 0.25T, k = 1,2, 3..., (see, Fig. [5] (a)). In this 
way, for the duration of pulses T = 0.57 -1 , and the time 
interval between them r = 47 -1 the Wigner function 
is depicted in Fig. IHb) for time measurement intervals 
t = kr + 0.8T. One might expect that this result would 
be approximately close to the pure |2) Fock state Wigner 
function displayed for comparison in Fig. [Hla). This 
Wigner function displays ring signature with the center 
at x = y = in phase-space. Indeed, in the general 
these results are qualitatively similar involving also neg- 
ative part, however, the cyclic symmetry of pure |2) has 
not displayed in state of PDAO that acquires two-fold 
symmetry in phase-space due to parametric excitation. 



V. CONCLUSION 

We have demonstrated production of various quantum 
states for the parametrically driven anharmonic oscillator 
in the regimes of low excitation and in complete consid- 
eration of dissipative effects. This investigation contin- 
ues our previous analysis [6j devoted to creation of Fock 
states as well as superpositions of Fock states in the spe- 
cific regime of periodically pulsed anharmonic oscillator 
for time-intervals exceeding the characteristic decoher- 
ence time. In this paper, we have proposed the other 
nonlinear oscillatory system that is excited by the train 
of Gaussian laser pulses through the degenerate down- 
conversion process. It can be realized as OPO combined 
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FIG. 4. The excitation number distributions for definite tim 
t = It- 1.8T (b); t = 2r - 0.4T (c); t = 2r (d); t = 2r + 0.6T 
the same time- intervals. The black circles indicate the negativ 
\h = 5, EQ,h = 10, r = 4 7 ~\ T = 0.57" 1 . 




FIG. 5. Time evolution of averaged excitation number (a) and 
populations of Fock states (b). The parameters are: A/7 = 
-10, x/l = 5 ' En h = 10-3, r = 4 7 -\ T = O.57- 1 

with phase modulation. Preparation of quantum states 
in PDAO has been stipulated by strong Kerr nonlinear- 
ity as well as by two-quanta resonant condition. We 
have studied the role of phase-localizing processes on pro- 
duction of oscillatory nonclassical states on the level of 
few excitation numbers. In this way, the production of 
the states approximately close to the superposition state 
I*) = 1 (|0) + 1 2}) and |2) Fock state have been de- 



3-intervals (gray points noted in Figs. 3(a)): t — 2r — 2T (a); 
(e). The Wigner functions corresponding these distributions for 
i parts of the Wigner functions. The parameters are: A/7 = —2, 



scribed. 



(a) (b) 




FIG. 6. The Wigner function for pure |2) state (a). The 
Wigner function for t = 2r — 0.8T (b), for the parameters: 
A/7 = -10, xh = 5, E ^h = 10.3, T = O.57" 1 , r = 47" 1 . 
The ranges of negativity are indicated in the black. 
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